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ABSTRACT 

QSOs are believed to be powered by accretion onto massive black holes (BHs). In this paper, 
assuming that each central BH in nearby galaxies has experienced the QSO phase and ignoring 
BH mergers, we establish a relation between the QSO luminosity function (LF) and the local BH 
mass function (MF). The QSOLF is jointly controlled by the luminosity evolution of individual 
QSOs and the triggering history of the accretion onto seed BHs. By comparing the time integral 
of the QSOLF with that inferred from local BHs, we separate the effect of the luminosity evolution 
of individual QSOs from the effect of the triggering history. Assuming that the nuclear luminosity 
evolution includes two phases (first increasing at the Eddington luminosity with growth of BHs 
and then declining) , we find that observations are generally consistent with the expected relation 
between the QSOLF and the local BHMF and obtain the following constraints on QSO models and 
BH growth: (i) The QSO mass-to-energy efficiency e should be > 0.1. (ii) The lifetime (defined 
directly through the luminosity evolution of individual QSOs here) should be > 4 x 10 7 yr. The 
characteristic declining timescale in the second phase should be significantly shorter than the 
Salpeter timescale rg p , and BH growth should not be dominated by the second phase, (iii) The 
ratio of obscured QSOs/AGNs to optically bright QSOs should be not larger than 7 at Mb ~ —23 
and 3 at Mb ~ —26 if e = 0.31, and not larger than 1 at Mb ~ —23 and negligible at Mb ~ —26 
if e = 0.1. (iv) It is unlikely that most QSOs are accreting at super-Eddington luminosities. We 
point out that it is hard to accurately estimate the value of the QSO lifetime estimated from the 
QSOLF and/or the local BHMF, if it is longer than a certain value (e.g., ~ 4rs p in this study). 
We discuss the importance of accurate measurements of the intrinsic scatter in the BH mass and 
velocity dispersion relation of local galaxies and the scatter in the bolometric correction of QSOs. 
We also discuss some possible applications of the work in this paper, such as to the study of the 
demography of QSOs and the demography of normal galaxies at intermediate redshift. 

Subject headings: black hole physics - galaxies: active - galaxies: evolution - galaxies: luminosity 
function, mass function - galaxies: nuclei - quasars: general 



1. Introduction 

The exploration of the relation of QSOs with massive black holes (BHs) in nearby galaxies has been of 
considerable interests since the discovery of QSOs. On the one hand, QSOs are suggested to be powered 
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by gas accretion onto massive BHs (Salpeter 1964; Zel'dovich & Novikov 1964). This (currently widely 
accepted) model suggests that a population of massive BHs as "dead" QSOs exist in nearby galactic centers 
(Lyndcn-Bell 1969). Furthermore, the total mass density of these remnant BHs and the typical BH mass 
in nearby galaxies can be inferred from the total energy density radiated in photons by QSOs (Soltan 1982; 
see also Rees 1984). These simple and elegant arguments motivate the search for BHs in nearby galaxies 
(e.g., Kormendy & Richstonc 1995; Magorrian et al. 1998; Gebhardt et al. 2003; Kormcndy 2003). On 
the other hand, as a result of the endeavor in the past two decades, not only has the existence of the 
massive dark objects (which are presumably BHs here) in most nearby galactic centers been confirmed, but 
also dramatic progress on their demography has been recently achieved (e.g., Kormendy & Gebhardt 2001; 
Tremaine et al. 2002 and references therein). Comparisons between the properties of local BHs and those 
inferred from the QSO model may shed new light on our understanding of the BH growth, the accretion 
physics, the mechanisms to trigger and quench nuclear activities, the formation and evolution of galaxies 
etc. (e.g. Cavaliere & Padovani 1989; Small & Blandford 1992; Fabian & Iwasawa 1999; Salucci et al. 1999; 
Elvis, Risaliti & Zamorani 2002; Marconi & Salvati 2002; Yu & Tremaine 2002; Fabian 2003). 

The comparison between the local BH mass density and the local energy density in QSO photons has 
suggested that local BHs acquire most of their mass through accretion during the QSO/AGN phase (Yu & 
Tremaine 2002; Aller & Richstone 2002; Fabian 2003). The mass distribution of local BHs as the remnant 
of nuclear activities is therefore controlled by the triggering history of the accretion onto seed BHs and 
the luminosity evolution of individual triggered nuclei with a mass-to-energy conversion efficiency e (i.e., 
the growth of individual BHs due to accretion). The triggering rate is usually believed to be related to 
the formation and evolution of galaxies and is a function of cosmic time, and the luminosity evolution of 
individual nuclei is believed to contain information on the accretion process in the vicinity of the BH and 
is a function of the physical time spent since the triggering of the accretion onto seed BHs. The triggering 
history and the luminosity evolution jointly control the QSO luminosity function (LF) as a function of 
luminosity and redshift; however, each cannot be uniquely determined from the QSOLF itself, partly as a 
result of the mixing of their effects on the QSOLF. In this paper, by using the local BH mass function (MF) 
as an additional constraint on QSO models and establishing a new relation between the QSOLF and the 
local BHMF, the luminosity evolution of individual QSOs is isolated from the triggering history of the QSO 
population in their effects on the QSOLF, and then we use observations to test QSO models and provide 
constraints on the QSO luminosity evolution and BH growth. 

A relation between the QSOLF and the local BHMF has been achieved by Yu & Tremaine (2002) under 
the assumptions that the seed BH mass is negligible and the luminosities of QSOs are only an increasing 
function of their central BH mass. In the study of this paper, we relax these two assumptions. In addition, 
the relation in Yu & Tremaine (2002) includes the effect of BH mergers. However, BH mergers are ignored 
in the relation established in this paper for the following reasons: (i) BH mergers are not shown to play a 
significant role or not necessarily required at least for growth of high-mass (> 10 8 M Q ) BHs, if e ~ 0.1 (Yu 
& Tremaine 2002); (ii) currently, the BH merger process and rate are very uncertain; and (iii) comparison of 
observations with the expectation obtained by ignoring BH mergers may also provide considerable insights 
in the role of BH mergers. The detailed difference between these two relations will be further discussed in 
this paper (see § 2.3 and 4.2.1). 

This paper is organized as follows. In § 2, by studying the continuity equation for the BH mass 
and nuclear luminosity distribution, we establish the relation between the QSOLF and the local BHMF. 
The triggering history of the accretion onto seed BHs is (implicitly) considered in the continuity equation 
but circumvented in the relation between the QSOLF and local BHMF. Only the luminosity evolution of 
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individual QSOs is incorporated in this relation. In § 3 we obtain the local BHMF by using the Sloan 
Digital Sky Survey (SDSS) observation results on both early-type and late-type galaxies and the empirical 
relations on the demography of galaxies and BHs (the BH mass and velocity dispersion relation, the Tully- 
Fisher relation, etc.). We also review the QSOLF obtained from large optical surveys. In § 4 we combine 
the observations with the relation obtained in § 2 to provide accurate observational constraints on the 
luminosity evolution of individual QSOs and BH growth, such as the QSO lifetime, the mass-to-energy 
conversion efficiency, the role of obscuration in BH growth, etc. The results are discussed in § 5. In § 6 
we further discuss some possible applications of the framework established in § 2, such as to the study of 
the demography of QSOs and the demography of normal galaxies at intermediate rcdshift. Finally, our 
conclusions are summarized in § 7. 

In this paper we set the Hubble constant as Hq = 100/ikms -1 Mpc -1 , and if not otherwise specified, 
the cosmological model used is (Qm, ^a, h) — (0.3, 0.7, 0.65) (Wang et al. 2000). 



2. The expected relation between the QSO luminosity function and the local BH mass 

function 

2.1. The continuity equation 

In this subsection we describe the evolution of the BH distribution by a continuity equation, which will be 
used to establish the relation between the QSOLF and the local BHMF. 

We define AT(ti, Mbh,0i L, t) (t > ti) so that jV(tj, Mbh,0: L, t)dtidM BBO dLdt is the comoving number 
density of local BHs with such properties: the nuclear activity due to accretion onto their seed BHs was 
triggered during cosmic time U — > tj + dtj, their nuclear luminosities were in the range L — > L + dL at cosmic 
time t, and these BHs have mass in the range Mbh.o - * Mbh,o + dMsH,o at present time to- We assume 
that the change rate of the nuclear luminosity L is a function only of (ti, Mbh.o, L, t). Ignoring BH mergers, 
we can use the following continuity equation to describe the evolution of Af(U, Mbh.o , L, t): 

dAf(U,M BHfi ,L,t) d[L(U,M B H,o,L, t)Af(U, M BH ,o, L, t)} 

dt dL ' U 

We define the nuclear LF n^(L,i) as follows: 

/>oo f-t 

n L (L,t)= dM BH ,o / dtiAf{ti,M BIlfi ,L,t) (2) 
Jo Jo 

so that n L (L,t)AL is the number density of the local BHs that had nuclear luminosities in the range L — ► 
L + dL at cosmic time t. 

Below we illustrate the relation of equation (1) with some continuity equations of the QSOLF or BHMF 
in the literature (e.g., Caditz & Petrosian 1990; Small & Blandford 1992). By integrating equation (1) over 
Mbh,o from to oo and over ti from to t, we may obtain the evolution of n^(L,t) as follows: 

dn L (L,t) , d[(L)n L (L,t)} _ a/T „ . . 

dt + dL ~ S{ - L ^ I 3 ) 

where (L)(L, t) is the mean change rate of L defined by 

. _ J^dMBH.o JodU L(ti,M BKfi ,L,tyf(ti,M B K,o,L,t) 
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and S(L, t) is the source function denned by 

/>oo 

S(L,t) = / M(ti,M BHfi ,L,t)\ ti=t dM BHfi (5) 
Jo 

and describing the triggering rate of nuclear activities of seed BHs. If L{t i: Mbh.Oj L, t) is a function only of 
L and t and riL(L,t) is replaced by the QSOLF, equation (3) will be identical to equation (1) in Cavaliere 
ct al. (1971) or equation (3) in Caditz & Petrosian (1990). If we replace L in equation (3) with BH mass of 
the progenitors of local BHs, equation (3) will look the same as equation (8) in Small & Blandford (1992) 
(which describes the evolution of the BHMF) . 

2.2. The time integrals of the nuclear/QSO luminosity function 

By integrating equation (1) over L from to oo, we have the conservation of the number density 

/•OO 

N(U,M BHfi ,t)= N(ti, Mbh,o> L, t)dL (6) 
Jo 



that is, 

dN(U,M B K,o,t) 
dt 

According to equation (7), we have 



0. (7) 



N(ti,M BBfi ,t) = JV(ti,M B H,o,to)., (8) 

For the same BH mass Mbh.o at present, we assume that their nuclear luminosity evolution is a function 
only of the age of their nuclear activities r = t — U. Thus, using equations (6) and (8), we have 

N(U,M BHfi ,L,t) = N(U,M BBfi ,t)S[L- £(Mbh,o,t)] 

= N{ti,M BHfi ,t )S[L-C{M BH ,o,T)}, (9) 

where S(x) is the Dirac function and C(M bb _q,t) is the nuclear luminosity of the progenitor of the BH 
A^bh.o at age r. As seen from equations (2) and (9), the nuclear luminosity evolution (incorporated in the 
5-function) and the nuclear activity triggering history [N(ti, Mbh,o, to)] jointly contribute to the nuclear LF. 
By first integrating equation (9) over tj from to t and over t from to to and then changing the integration 
variables (t,U) to (U,t), we have 

f dt f dtiN{t u M BBt o,L,t) = [°dt[ dU N(U,M BBfi ,to)5[L-jC(M BIlfi ,T)] 
Jo Jo Jo Jo 

= f dU N(U,M BBfi ,to) f dr S[L- C(M bb , ,t)] 
Jo Jo 

= f dti N(ti, M B h,0j to) TTtTm 1 ' 

Jo V |d£(M B H,0,T)/dT| r=Tfe | 



(10) 



where Tk(L, Mbh,o) {k — 1, 2, ...) are the roots of the equation £(Mbh,o, t) — L = (0 < t < to — ti). 1 We 
assume that the nuclear activities of all the local galaxies are quenched at present. Thus, for local BHs with 



1 In equation (10), it is assumed that d£(A4"BH,0i r)/dr ^ at r = rj, (k = 1,2, ..). Formulae (10), (11), (14) and (15) are 
not difficult to generalize even if d£(MBH,o> t)/At = at r = r*,, and other formulae will not be changed. 
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mass Mbh.o, they have the same roots Tk(L, Mbh.o) {k — 1, 2, ...) for the same L since they have experienced 
the same evolution of £(Mbh,Oj'?~) before their quenching, even though their accretion onto their seed BHs 
may be triggered at different time U. Hence, the sum term J2k m equation (10) does not depend on U, and 
we have 

/•to rt i 

I dt dtj Af(U, M B h,o, L, t) — um bh (M BH ,o,to)> tt^ttj rjT-, r, (H) 

Jo Jo \dL{M BHfi ,T)/dT\ T = Tk \ 

/•to 

n MBH (M B K.o,to) = / N(t u M BB .o, t )dU (12) 
Jo 

is the local BHMF and nM BH (-^bh,o, io)dM B H,o gives the number density of local BHs with mass in the 
range Mbh,o — * Mbh,o + dMsH,o- The lifetime of the nuclear activity for the BH with current mass Mbh,o 
can be expressed by 



where 



7iifo(M B H,o) = / dr (13) 



/ 



dL E| d£( M B H,o,r)/dT| r=r J- (U) 



Note that the integration in equation (13) is over the period when the nucleus is active, and equation (14) is 
also restricted to the luminosity £(Mbh,o, t ) that is taken as active (e.g., higher than a certain luminosity 
limit). The definition of "active" of galactic nuclei may be different in different contexts, and equations 
(13) and (14) (such as their integration limits) may be adjusted to appropriate forms according to different 
definitions. During the period of the nuclear activity, the fraction of the time (or the probability) with 
luminosity in the range L — > L + dL can be given by 

P(L|MBH <° )di S ^fihd ? |d£(M B H,,V|_J (15) 

with 

J P(L\M Bii ,o)dL = 1. (16) 
Applying equations (14) and (15) in equation (11), we have 

/■to rt 

dt dt i ^(t i ,M BKfi ,L,t) = nM B AM BHi o,to)ni{ e (M B ii,o)P{L\M BHi0 ). (17) 
Jo Jo 

By integrating equation (17) over Mbh,o and using equation (2), we have the time integral of the nuclear 
LF as follows: 

pt() poo 

/ n L (L,t)dt= dM BH ,o n MBH (MBH,o,*o)'nifo(MBH,o)i 5 (i|MBH,o) = T L , local (L, t )- (18) 
Jo Jo 

We define the QSOLF ^>L(L,t) so that ^l(L, t)dL is the comoving number density of QSOs with 
luminosity in the range L — > L + dL at cosmic time t, and we define the time integral of the QSOLF as 
follows: t 

Tl.qsoOMo) = [ * L (L,t)dt. (19) 
Jo 
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The QSOLF in equation (19) includes the contribution from both optically bright QSOs and any other QSOs 
that are obscured or not seen in optical bands but might be detectable in other bands (see Fabian 2003). 
In this paper ^L{L,t) denotes the LF of the "live" QSOs in the distant universe and m,(L,i) denotes the 
nuclear LF of the progenitors of the local BHs. Based on the cosmological principle, similar to Soltan's 
argument (1982), the QSOLF ^l(L, t) represents the evolution of the nuclear luminosity of local BHs, i.e., 

n L (L,t) = V L (L,t). (20) 

Thus, with equations (18)-(20), we have 

7l,QSC>(£, to) = 7t,local(i, to). (21) 

The physical meaning of equations (18) and (21) can be understood clearly as follows. For the progenitor 
of each local BH with mass Mbh,o, the average time that it has spent in the nuclear luminosity range 
L — > L + dL is T\n c (MBH y o)P(L\MBH,o) , an d taking QSOs as the progenitors of the local BHs, the total time 
spent in the range L — > L + dL by the progenitors of all the local BHs with mass Mbh.o i n a unit comoving 
volume should just be the time integral of ^l(L, t)dL. We note that the quantity of the time integral of the 
QSOLF has been (implicitly) used before in the literature, such as in some BH mass density relations between 
local BHs and QSOs or the definition of the QSO mean lifetime, which will be further discussed in § 2.3, 
and we also note that Blandford (2003) points out that a simple model of accretion implies a quantitative 
relationship between the time integral of the QSOLF and the local BHMF. Relation (21) established above 
on the time integral of the QSOLF will be the base to constrain the luminosity evolution of QSOs and BH 
growth in this paper. 

In addition, note that in equation (21) or (18) the luminosity L is assumed to be only a function of 
-Mbh,o and r. For the more complicated case that L also depends on some other parameters, we may take the 
lifetime and probability function in equation (18) as the averaged result over other parameters or generalize 
these equations by including other parameters. 



2.3. Mean lifetime, and total/partial BH mass densities 

Below we show with appropriate assumptions, how the definition of the mean QSO lifetime and some relations 
on the BH mass density obtained in the literature (e.g., Soltan 1982; Yu & Trcmaine 2002; Haiman, Ciotti 
& Ostriker 2003) can be incorporated in the framework established above: 

• Mean QSO lifetime: the lifetime riif (MBH.o) in equation (18) also represents the lifetime of "live" 
QSOs whose central BH masses will be Mbh,o at present. By integrating equation (21) over L and 
using equation (16), we have 

//■to roc 
dL / V L (L,t)dt = / dM BH ,0 ^M BH (-^BH,0)to)7"life(-^BH,o)- (22) 
Jo Jo 

We may define the mean lifetime of QSOs as follows: 

= / °° dM BH ,on MBH (Mbh.o, *o)7iife(M B H,o) = I dL Jq° ^L(L,t)dt 

/ 00 dMBH,o"MBH(^BH,o,to)ff[Tiifo(M B H.o)] ~~ J °° dM BHi „«M B H (M B n,o, to)fl"hife(M BH ,o)] ' 

(23) 

where equation (22) is used and H{x) is the Heaviside step function defined by H(x) = 1 if x > and 
H(x) = otherwise. We denote the mass of the progenitors of the local BHs when they had a nuclear 
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luminosity L (or the central BH mass of a QSO with luminosity L) as Mbh- If L is an increasing 
function only of Mbh and only the period with luminosities brighter than a certain value L(M BH ) is 
taken as active phase, we have riif e (MBH,o) > if Mbh,o > M BH and riif e (MBH,o) = otherwise. Thus, 
the mean lifetime of QSOs is given by 

J m / jh ^Mbh(Mbh,o, to) dM BH ,o 

which is identical to equation (59) in Yu & Tremainc (2002; see also cq. [24] in Haiman, Ciotti & 
Ostriker 2003). 

Total BH mass densities: we denote the QSO bolometric luminosity produced by a mass accretion rate 
M acc as Lboi = eM acc c 2 = cMbh/(1 — e ), where Mbh = M acc (l — e) is the growth rate of BH mass 
and c is the speed of light. (The subscript "bol" represents the bolometric luminosity, and the symbol 
L in this paper may be either the bolometric luminosity or the luminosity in a specific band, if not 
specified.) By multiplying equation (21) by (1 — e)L^ \/ (ec 2 ) = Mbh and then integrating over L from 
a given value L' to oo, we have 

roo /i \j rta poo poo 

/ dL y U^L * L (L,t)dt= dM BH ,o n MBH (Mbh.o, to) / dL M BH r li{e (M B H,o)P(L\M BH . ). 

J L> ec JO Jo J L' 

(25) 

By setting L' = 0, equation (25) becomes 

f°° (1 — e)Lh f to r°° 

/ dL ± >_*5± / * L (L,t)dt= dM BH ,o M B H,on MBB (M BKtQ ,to), (26) 

Jo ec Jo Jo 

where 

/ dL MBHTlife(M B H,o)^(i|MBH,o) - / M BH dt = M BH .o - (M BH ,i> (M BH) o) (27) 
Jo Jo 

is used; (Mbh, i) (Mbh, o) is the average mass of the seed BHs that will have mass Mbh.o at present and 

is ignored in equation (26). Equation (26) is identical to Soltan's argument (1982) relating the total 

local energy density in QSOs to the total BH mass density in nearby galaxies. 

Partial BH mass densities: if the QSO luminosity L is only an increasing function of the BH mass 
M B h 5 we have 

/oc 
dL M BH Tli f c(MBH,o)P(i|MBH,o) 

Mbh.o - max[M BH (L'), M3H, 4 (M BH ,o)], if M BH ,o > max[M BH (i'), Mbh, i (Mbh, o)] (2g) 
0, otherwise 

where the seed BH mass Mbh,i is assumed to be a function only of Mbh,o- By applying equation (28) 
in equation (25) and assuming Mbh, i (Mbh, o) < M BB (L') 7 we have 

f d£ (l-e)£boi f\ L[L t)At= f dMBH f) [Mbh,o-Mbh(£>m bh (Mbh,o,*o). (29) 

J L> ec Jo JMbh(L') 

Actually, even if the seed BH mass is not a function only of Mbh,o, equation (29) will still hold at 
a given L' as long as all the seed BH masses of the local BHs M B h,o[> M bb (L')] are smaller than 
M BB {L'). Equation (29) is identical to the relation between the partial BH mass density in nearby 
galaxies and that accreted during bright QSO phases in Yu & Tremaine (2002) if BH mergers are 
ignored (see eq. 30 in Yu & Tremaine 2002). 
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2.4. £(M bh ,o,t) 

Given the QSOLF and local BHMF, the nuclear luminosity evolution [or Tiif e (MBH,o)-F > (-k|-^BH,o)] cannot 
be uniquely determined from equations (18)-(21). Below we use some physical arguments to assume a form 
of the luminosity evolution £(Mbh\o> t)- The parameters involved in the assumed form will be constrained 
by observations in § 4.2. 

It is generally believed that the growth of a BH involves two phases after the nuclear activity is triggered 
on (see the discussion on the "feast and famine" model in Small & Blandford 1992 and Blandford 2003). 
In the first (or "demand limited") phase, there is plenty of material to supply for the BH growth; however, 
not all of the available material can contribute to the BH growth at once and the BH growth is limited by 
the Eddington luminosity. With the decline of material supply, the BH growth enters into the second (or 
"supply limited" ) phase and the nuclear luminosity is expected to decline below the Eddington luminosity. 
For simplicity, we assume that the two phases appear only once in this paper, although this assumption 
is not required in relation (21). The possibility of more complicated accretion patterns deserves further 
investigation by applying corresponding models in relation (21). 

In the first phase, after the BH is triggered at time t i7 we assume that it accretes with the Eddington 
luminosity for a time t\ and its mass increases to be Mbh,i at time t = U + t\ = t\. The mass-to-energy 
conversion efficiency e is assumed to be a constant. Thus, the nuclear luminosity in the first phase increases 
with time as follows: 

£boi(r) - LEdd(M BH .i)exp ( — — ) , 0<r<n, (30) 

V T s P / 

where L~Edd{MBn) is the Eddington luminosity of a BH with mass Mbh and 



r Sp = 4.5 x 10 7 



0.1(1-6) 



yr (31) 



is the Salpeter time (the time for a BH radiating at the Eddington luminosity to e-fold in mass). In the 
second phase, we assume that the evolution of the nuclear luminosity declines as follows (e.g., Hachnelt, 
Natarajan & Rees 1998; Haiman & Loeb 1998): 

£ bol (r) = { L Edd(M BH ,i) exp , for n < r < n + £r D , (32) 

[ 0, forr>ri+^r D , 

where td is the characteristic declining timescale of the nuclear luminosity. We assume that QSOs become 
quiescent when the nuclear luminosity declines by a factor of 77 = cxp(— £) compared to the peak luminosity 
-^Edd(-^BH,i), so there is a cutoff of the nuclear luminosity at r = t\ + £td in equation (32). The factor 
£ is set to — ln(10 -3 ) = 6.9 here, since after decreasing by a factor of ry = 10 , the nuclear luminosity of 
BHs even with a high mass <~ 10 9 M Q will become fainter than the luminosity range (Mb < —20 in § 4) of 
interest in this paper. With the assumption that all QSOs are quenched at present (i.e., i — t.i — t\ ^ td), 
the BH mass at present is given by 

M BH; o - M BH ,i + f {1 ~ € l Lhol dt =(1+™-) Mbh.l (33) 
J tl « 2 V t Sp J 

In equation (33), the efficiency e is assumed to be the same as the efficiency in the first phase and not to 
change with the decline of the nuclear luminosity. It is important to know how the realistic efficiency evolves 
with the change of the nuclear luminosity or other parameters (which would depend on the evolution of both 
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the accretion rates and BH spins; a detailed study about this is beyond the scope of this paper). If the 
nuclear luminosity of a BH is smaller than its Eddington luminosity by a factor of 10 or more, according to 
current accretion models, the BH might accrete material via the advection-dominated accretion flow (ADAF) 
with low efficiency e <C 0.1 (e.g., Narayan & Yi 1994), rather than via the thin-disk accretion with efficiency 
e <~ 0.1 — 0.3 near the Eddington luminosity. Considering this possibility of the very low efficiency at the late 
stage of the luminosity evolution, the contribution from the very low efficiency stage to the time integral of 
the nuclear LF (in the luminosity range Mb < — 20 shown in this paper) is insignificant and our conclusions 
will still hold for other stages of the nuclear activity (which can be inferred from the results obtained in 
§ 4.2). We also assume that e is irrelevant to the BH mass Mbh,o in this paper. 

By using equations (30)-(33) in equations (14) and (15), we have the QSO lifetime given by 

Tlife = n + £t d , (34) 

and the BH with present mass M B h,o has such a probability distribution of the nuclear bolometric luminosity 
in its evolution history: 

P(£boi|M BH ,o) - — (fi^ + fo-^-) , (35) 

Tlife \ ^bol ^bol/ 

^ _ f 1, if iEdd(M B H,i)exp(-Ti/T S p) < Lboi < L Edd (M B n,i) ,^q\ 
\ 0, otherwise 

/d = f 1, if L Edd (M BH j) exp(-£) < L bol < £ E dd(M BH ,i) ^ 
1 0, otherwise 



where 



and 



In the analysis above, the timescales t\ and td are not necessarily constants, and they may be a function 
of Mbh,o- The dependence of t\ on Mbh,o would depend on the distribution of seed BHs, which are poorly 
known and beyond the scope of this paper. Below we always assume that t\ and t b are irrelevant to Mbh,o- 



3. The local BH mass function and QSO luminosity function obtained from observations 

3.1. The local BH mass function um bh (Mbh,o, h) 

In this subsection we use the velocity-dispersion distribution of local galaxies (§ 3.1.1) and the empirical BH 
mass-velocity dispersion relation (§ 3.1.2) to obtain the local BHMF (§ 3.1.3). 



3.1.1. Local velocity-dispersion function n a (a,t ) 

We define n a (a,t ) as the velocity-dispersion function of the hot stellar components of local galaxies so 
that n rT (a,t )da represents the comoving number density of local galaxies in the range a — > a + da (by 
"hot" component we mean either an elliptical galaxy or the bulge of a spiral or SO galaxy). The velocity- 
dispersion distribution n a (a,to) includes the contribution by both early- type galaxies n® arly (cr, to) and late- 
type galaxies n^ te (a, to), that is, 

n„(*,to) = < arly Mo) +4 ate (a,t ). (38) 
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Recent study on a sample of nearly 9000 nearby early-type galaxies obtained by the SDSS has provided 
the velocity-dispersion distribution in early-type galaxies as follows (see eq. 4 in Sheth et al. 2003, and 
Bcrnardi ct al. 2003): 

early, . x , ( <? V [~W^] (3 ,„ Q . 

n a y {oM) = fa\ — \ , (39) 

Vcr*/ r(a//i) cr 

where the best-fit values of (</>*, cr„, a, /3) are (0.0016 ± 0.0001, 88.8 ± 17.7, 6.5 ± 1.0, 1.93 ± 0.22), fa is the 
comoving number density of local early-type galaxies in units of Mpc~ 3 , and cr* is in units of kms -1 . 

We obtain the velocity-dispersion distribution in local late- type galaxies njf te (cr, to) in the following 
steps: (i) Following Sheth et al. (2003), we get the LF of the late-type galaxies by subtracting the LF of the 
early-type galaxies (Blanton et al. 2003) from the LF of total galaxies (Bcrnardi et al. 2003). (ii) Following 
Sheth et al. (2003), we get the distribution of the circular velocity v c in late-type galaxies by using the LF 
of the late-type galaxies obtained above and the following Tully-Fisher relation (Giovanelli et al. 1997): 

log(2w c ) = 1.10 - (Mi - 51og/i)/7.95, (40) 

where M t is the absolute magnitude of the galaxies in the / band, accounting for the intrinsic scatters around 
relation (40) and the inclination effects of galaxies (see details in Sheth ct al. 2003). (iii) We get the velocity- 
dispersion function of late-type galaxies by using the circular-velocity distribution of the late-type galaxies 
obtained above and the following correlation between the circular velocity and the velocity dispersion of the 
bulge component (see eq. 3 in Baes et al. 2003; see also Ferrarese 2002): 

l0g (200^) = ( °- 96 ± °- n > bg (200i^i=l) + ( °- 21 ± °-° 23) - (41) 

The intrinsic scatter of the correlation given by (41) is small (< 0.15 dex; see Fig. 1 in Baes ct al. 2003) and 
ignored in our calculation. Relation (41) may not hold for a < 80kms _1 , which corresponds to BH mass 
< 4 x 10 6 M according to the BH mass- velocity dispersion relation below (cq. 42) and is beyond the range 
that we focus on in § 4.2. 



3.1.2. The BH mass and velocity dispersion relation 

Studies of central BHs in nearby galaxies have revealed that the BH mass and the velocity dispersion of the 
hot stellar component of the host galaxy follow a tight correlation (Trcmaine et al. 2002; Ferrarese & Merritt 
2000; Gebhardt et al. 2000). The logarithmic of the BH mass at a given velocity dispersion o has a mean 
value given by (Tremaine et al. 2002) 

(log(M BH ,o)|crc> = A + 7 log( C r c /200kms- 1 ), (42) 

where Mbh\o is in units of M , 7 = 4.02 ± 0.32, A = 8.18 ± 0.06 have been adjusted to our assumed Hubble 
constant h — 0.65 (see section 2.2 in Yu & Tremaine 2002), and er c is the luminosity- weighted line-of-sight 
velocity dispersion within a slit extending to the effective radius. Note that a c in the correlation (42) is the 
velocity dispersion within a slit extending to the effective radius R (Gebhardt et al. 2000), while a in the 
SDSS (eq. 39) is the velocity dispersion within a circular aperture extending to R a /8. However, replacing a e 
with a in equation (42) will not cause much difference as the two definitions should give very similar results 
(Trcmaine et al. 2002). Merritt & Ferrarese (2001) give an alternative version of the correlation (42) with a 
steeper slope, 7 = 4.72. The reasons for this difference in slopes are discussed by Tremaine et al. (2002). We 
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perform the calculations using both versions of the correlation and the difference on the local BHMF will be 
discussed in § 3.1.3. 

Note that relation (42) is fitted in the logMBH,o~log<7 space. We assume that the distribution in 
log Mbh, o at a given a is Gaussian with intrinsic standard deviation Ai og M B H,o; which is independent of a 
and thus can be written as 

(logMBH.O- (lQgM B H.o|g)) 2 l 

log M BH ,o 

According to Trcmaine et al. (2002), the intrinsic scatter in log Mbh, o should not be larger than 0.25-0.3 dex. 
Nevertheless, to check the effect of Ai og M BHi0 on the results, below we show the results obtained with a value 
of A, oe 

m bh ,o(~ 0.4 dex) higher than 0.27 dex, as well as those obtained with Ai og M BH () = and 0.27 dex. 

With the velocity-dispersion function and the BH mass-velocity dispersion relation (cqs. 42 and 43), 
the local BHMF is given by 

t ) = In 10 [ A^BH,o-P(logM B H,o|o')n (T (CT, t )d<T. (44) 



P(logM BH ,o|CT) = 



2ttAi, 



■ exp 



OgM B H,0 



3.1.3. Results: nM BH (^bh,o, *o) 

Using the velocity-dispersion function in § 3.1.1, the BH mass-velocity dispersion relation in § 3.1.2, and 
equation (44), we obtain the local BHMF riM Bii (Mbh.o, ^o) an d show the distribution of Mbh.o^ a/bh (^bHi ^o) 
in Figure 1. In Figure la, different solid lines represent the results obtained with different Ai G g Mbh.o (=0, 
0.27, and 0.4 dex from bottom to top at the high-mass end, respectively, see eq. 43). As seen from panel 
(a), the local BHMF at the high-mass end (> 3 x 10 8 M Q ) is significantly affected by the intrinsic scatter 
A} g Mbh o • For example, the BHMF obtained with A} g Mbh o = 0.27 dex (middle solid line) is larger than 
that obtained with Ai og M B H,o = (bottom solid line) by a factor of 5 at Mbh.o — 1O 9 M and by a factor 
of more than 10 at Mbh,o — 4 x 10 9 M q and if Ai og M BH>[) = 0.4 dex, the BHMF (top solid line) is larger 
than that obtained with A] og M BH = by a factor of 8 at Mbh,o — 10 9 M Q and a factor of more than 100 
at Mbh.o — 4 x 10 9 M Q . The solid lines are the results obtained by setting the variables (</>*, <r*, a, /?) in 
equation (39) to the mean of the best-fit values. We also show the uncertainty of the local BHMF due to the 
1-cr error of the fitting values in the velocity-dispersion function of early- type galaxies in panel (a). The two 
dotted lines adjacent to each solid line represent the result after considering the 1-a error of the mean value of 
(cr*,a,/3) in equation (39), which is obtained by setting (cr„, a, /3) = (88.8, 6.5, 1.93) ± (+17.7, -1.0, +0.22). 
Note that the signs set to the errors of the variables are not all the same because the best-fit values of 
(<r*,a,f3) are strongly correlated with one another (see fig. 4 in Sheth et al. 2003). We do not consider the 
error in <f>„ here since it is small and the uncertainty caused by it is negligible compared to that caused 
by the error of other variables. As seen from panel (a), the uncertainty of the local BHMF due to the 
uncertainty in the velocity-dispersion function of early-type galaxies is negligible compared to that due to 
the intrinsic scatter Ai og M B Ho m the Mbh.o — 17 relation. In panel (a), the dashed lines show the local 
BHMF in late-type galaxies also obtained with A] g Mbh o —0, 0.27 and 0.4 dex from bottom to top at the 
high-mass end, respectively. So far, it is hard to give an accurate estimate of the uncertainty in the estimate 
of the velocity-dispersion function of late-type galaxies here. We believe that the effect on the local BHMF 
due to the uncertainty in the velocity-dispersion function of the late-type galaxies can be ignored at the 
high-mass end (>4x 10 7 M ), where the local BHMF is dominated by early-type galaxies. The local BHMF 
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is dominated by late-type galaxies at the low-mass end (< 10 7 M Q ). We expect that the local BHMF with 
BH mass in the range from 4 x 1O 6 M0 to 4 x 1O 7 M has considerable accuracy, for example, within 50%. 

In Figure lb we show the effect on the local BHMF due to different versions of the slope 7 in the 
Mbh,o — & relation (eq. 42). The solid lines are the same as those in panel (a), which are obtained with 
7 = 4.02 (Tremaine et al. 2002). The dot-dashed lines are obtained with 7 = 4.72 (Merritt & Ferrarese 
2001). As seen from panel (b), the difference of the BHMF at the high-mass end due to the difference in 7 
is not so significant as that due to the difference of the intrinsic scatter Ai og jf BH in the Mbh,o — a relation. 
Below we always show the results obtained with 7 = 4.02 (Tremaine et al. 2002). 

In principal, the local BHMF can also be obtained by using the luminosity distribution of the hot 
stellar components of local galaxies and its correlation with central BHs (e.g., Kormendy & Gebhardt 2001). 
However, we do not do so in this paper, because the total local BH mass density obtained by using the BH 
mass - luminosity relation appears to be larger than that obtained by using the BH mass - velocity dispersion 
relation by a factor of more than 2, and the reason for this difference is not yet clear (Yu & Tremaine 2002). 
As argued in Yu & Tremaine (2002), we believe that the results obtained by using the BH mass - velocity 
dispersion relation are more reliable, as the correlation between the BH mass and the velocity dispersion is 
tighter and the result does not depend on the uncertain bulge-disk decomposition. 

Note that although the local BHMF at the high-mass end is significantly affected by the intrinsic scatter 
in the Mbh.o — a relation, the total local BH mass density is not, which is higher than the total density 
obtained by setting Ai og M B H,o = only by a factor of exp[i(Ai og M BH ,o hil0) 2 ](=1.2 if Ai og M BH . = 0.27dex 
and 1.5 if Ai og M B H,o = 0.4dex (see cq. 12 in Yu & Tremaine 2002). 



3.2. The luminosity function of optically bright QSOs 

The LF of optically bright QSOs is often fitted with a double power law: 

^M B ( M B,z) = 1q0A{i3i + 1)[Mb _ m , b{z)] ™ w oaU3 2 + 1)[M b -M'(z)]> ( 45 ) 

where the superscript "opt" represents optically bright (or unobscured) QSOs, ^ (Mb, z)dMs is the 
comoving number density of QSOs with absolute magnitude in the range [Mb, Mb + dMs] at redshift z, 
and we have 

V Mb (M b , z) = * Lb (L b , z)\dL B /dM B \ = 0.92L B ^(L B , z). (46) 

Pei (1995) uses this form to fit the data set from Hartwick & Schade (1990) and Warren, Hewett, & Osmer 
(1995) on the basis of more than 1200 QSOs over the range of redshift 0.1 < z < 4.5. In the cosmological 
model of (Q m , Q A , ft) = (1, 0, 0.5), the QSOLF with absolute magnitudes -30 < M B < -23 can be fitted by 



equation (45) with the following parameters: 

* M = 6.7 x 10~ 6 ft 3 Mpc~ 3 mag" 1 , (47) 

M* B (z) = M* B (0) + 1.251og(l + z)- 2.5(fci2 + k 2 z 2 ), (48) 

M B (0) = -20.83 + 5 log ft, An = 1.39, k 2 = -0.25, (49) 

/3i = -1.64 and (3 2 = -3.52. (50) 



Note that the cosmological model in which the parameters above are fitted is different from the model 
adopted in this paper. In our calculations below, we have transferred the parameters in equations (47)-(50) 
into the cosmological model used in this paper, (fl m , ft) = (0.3,0.7,0.65). 
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M BH.O ( M ©) 



Fig. 1. — Local BHMF obtained from observations of the velocity dispersion distribution of nearby galaxies 
and the BH mass-velocity dispersion relation, (a) Solid lines represent the Mbh,o^m bh (-^bh,0)^o) m both 
early-type and late-type galaxies (sec cq. 38-44); dashed lines represent the Mbh.o"m bh {M-qk.o, ^o) m late- 
type galaxies. Different solid or dashed lines are obtained with different values of the intrinsic scatter 
AiogM BH ,o (=0, 0.27, and 0.4 dex from bottom to top at the high-mass end, respectively) in the Mbh,o — o 
relation (eq. 42). The two dotted lines adjacent to each solid line give the effect of the l-cr error of the 
best- fit parameters in the velocity-dispersion function of early- type galaxies (eq. 39). (a) shows that the 
local BHMF at the high-mass end (> 3 x 10 8 M Q ) is significantly affected by Aiog Mbh,o > an d the- uncertainty 
of the BHMF due to the uncertainty in the velocity-dispersion function of early-type galaxies is negligible 
compared to the uncertainty due to the intrinsic scatter Ai og M BH in the Mbh,o — v relation, (b) Solid lines 
are the same as those in (a), which are obtained with the slope 7 = 4.02 in relation (42) (Tremaine et al. 
2002), while the dot-dashed lines are obtained with 7 = 4.72 (Merritt & Ferrarese 2001). The difference of 
the BHMF at the high-mass end due to the difference in 7 is not as significant as that due to the difference 
of the intrinsic scatter Ai og M B H.o m the Mbh.o — f relation. See details in § 3.1.3. 
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Boyle et al. (2000) use the function form of equation (45) to fit a much larger data set from the 2dF 
QSO redshift survey (Boyle et al. 2000) and Large Bright QSO survey (Hewett, Foltz & Chaffee 1995), which 
contain over 6000 QSOs, and give the QSOLF with absolute magnitudes —26 < Mb < —23 and redshift 
0.35 < z < 2.3 in our standard cosmological model (£l m , Q\) = (0.3, 0.7) by the following parameters: 



V* M = 2.9 x lO^Mpc^mag- 1 , (51) 

M* B (z) = M B (0) - 2.5(M + k 2 z 2 ), (52) 

M B (0) = -21.14 + 5 log h, h = 1.36, k 2 = -0.27, (53) 

/3i = -1.58 and 2 = -3.41. (54) 



The quadratic dependence of the characteristic magnitude M B {z) on z in equations (48) and (52) shows 
an increasing characteristic luminosity with increasing redshift at low redshift (z < 2.5) and a decline of the 
characteristic luminosity at higher redshift, which is suggested by observations (e.g. Shaver et al. 1996), 
although the LF is not yet accurate enough to confirm the decline at z > 2.5. 

The QSOLF over the range 3.6 < z < 6 provided in Fan et al. (2001, 2003) gives a flatter bright-end 
slope [fii — —2.5) than equations (50) and (54); however, in our calculation below, we simply extrapolate 
equations (50) and (54) to high redshift because the detailed QSOLF at z > 3.5 docs not affect our results 
much (see Fig. 2 below; if we use the LF in Fan et al. 2001, our results change by less than a few percent). 

Using equations (45)-(54), we obtain the time integral of the LF of optically bright QSOs 

/•oo 

*]3b,qso( M *'*°)= / *T B , Q so( M B,t)dt = OML B T^ QSO (L B ,to) (55) 
Jo 

as a function of the absolute magnitude Mb shown in Figure 2. Both of the results obtained from the QSOLFs 
given by Boyle et al. (2000; shown as solid lines; eqs. 47-50) and those given by Pei (1995; shown as dotted 
lines; eqs. 51-54) are presented. In the calculation, the QSOLFs are extrapolated to the parameter space not 
covered by the observations. In Figure 2a we show the integral of the QSOLF over time < z < 2.3 as the 
bottom solid and dotted lines and the integral of the QSOLF over the whole cosmic time qso(-^B' *o) 
as the top solid and dotted lines. As seen from panel (a), the difference of the time integral of the QSOLF 
due to different versions of QSOLFs is negligible. In Figure 2b we show the fraction of 7"^* q SO (Mb, io) 
contributed by QSOs in the range < z < 2.3 (bottom solid and dotted lines), which is more than 50% 
at Mb > —26 and 75% at Mb > —23, and we also show the fraction contributed by QSOs in the range 
< z < 3.5 (top solid and dotted lines), which is more than 95% at any Mb in the figure. Thus, the effect 
of the uncertainty in the QSOLF at high redshift on the time integral of QSOLF can be neglected. 



4. Is the local BH mass function consistent with the QSO luminosity function? 

With the goal to obtain accurate observational constraints on the QSO luminosity evolution and the fun- 
damental parameters of the QSO model (see § 2.4), we compare the time-integral of the QSOLF with that 
predicted from local BHs in this section (see cq. 21 in § 2.2 and inequalities 60 and 62 below). Note that 
the luminosity in the assumed model on the QSO luminosity evolution in § 2.4 is the bolometric luminosity 
iboi- However, the QSOLF obtained from observations is usually in a specific band (which is the B band in 
this paper) . Before the comparison, we first transfer the bolometric luminosity to the luminosity in a specific 
band. 
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Fig. 2. — Time integral of the LF of optically bright QSOs as a function of the absolute magnitude Mb (see 
eq. 55). The results obtained by using the QSOLFs given by Pei (1995) (eqs. 47-50) are shown as dotted 
lines, and those obtained by using the QSOLFs given by Boyle et al. (2000) (eqs. 51-54) are shown as solid 
lines. The redshift range of the QSOLF is only 0.35 < z < 2.3 for Boyle ct al. (2000) and 0.1 < z < 4.5 for Pei 
(1995). In the calculations, both are extrapolated to the parameter space not covered by the observations, (a) 
Top solid and dotted lines represent the integral of the QSOLF over the whole cosmic time ^~Mb,QSo(^b, to)', 
bottom solid and dotted lines represent the integral of the QSOLF over time in the redshift range < z < 2.3. 
The difference due to different versions of the QSOLFs is negligible, (b) Top solid and dotted lines give the 
fraction of T^b QSo(-^ s ' *o) contributed by the QSOs with redshift < z < 3.5; bottom solid and dotted 
lines give the fraction contributed by the QSOs with redshift < z < 2.3. This figure shows that the effect 
due to the uncertainty of the QSOLF at high redshift on the time integral of QSOLF can be neglected. See 
§ 3.2. 
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4.1. The bolometric correction 

We denote Cb as the bolometric correction in the B band, defined through Lboi = CbL„ b , where L VB is the 
energy radiated at the central frequency of the B band per unit time and logarithmic interval of frequency. We 
assume that Cb is independent of the cosmic time t. We denote P{Cb l^boi) as the probability distribution 
function of the bolometric correction Cb at a given Lboi- Thus, we have the LF of the progenitors of local 
BHs in the B band as follows: 

n LvB {L VB ,t)= J C B P(C B \L hoi )n Lbo] (L hoh t) dC B - (56) 
With equations (18) and (56), we have 

Tl„ b ,\oca\{L VB , to) = J n L„ B {L VB , t) dt — J C_B-P(Cs |Lbol)TL bol Jocal(ibol, to) ACb- (57) 

Similar to equation (21), we have 

TL VB ,qso{L VB ,to) =T LuBt \ oca \{L VB ,to); (58) 
and according to equation (55), we have 

1~M B ,QSC>(MB,h) = TM B! iocal(M_B, to). (59) 

Considering that QSOs in the equation above include both optically bright QSOs and obscured QSOs, we 
have 

T Ms,QS0( M B, to) < T Mb , local (Mfl, to). (60) 

The ratios of the time integral of the QSOLF and that predicted from the local BHMF are given as follows: 

^B,< a ) S ^ 8 '*'! =l (61) 

1 M B , QSO{M b , to) 

and 

n^\M B ^)^ ^^p^ >l. (62) 

The value of lZ opt (MB, to) — 1 gives the number ratio of obscured QSOs to optically bright QSOs. 

Elvis et al. (1994) study the spectral energy density distribution of a sample of 47 QSOs and find that 
the bolometric corrections in the B band have a mean value Cb = 11-8 with standard deviation Ac B = 4.3. 
The origin of the non-uniform of the bolometric correction (or the spectral energy distribution) is not yet 
clear, which could be due to the difference in the accretion rate, the BH mass and spin, etc. For example, 
with increasing BH mass, the bolometric correction in the B band may decrease since the peak of its spectral 
energy distribution may move from a short wavelength toward the B band or long wavelength, and with 
increasing accretion rate, the bolometric correction in the B band may decrease since it may have a softer 
spectrum in the X-ray band. In this paper, for simplicity, we assume that the probability distribution 
function of the bolometric correction P(C_B|Lboi) is independent of Lboi and follows a Gaussian distribution 
as follows: 

(C B - C B ? 



P{C B \L hol ) = -=- cxp 



(63) 



In the study below we use two values of A Cb = 4.3 and [P (Cb\L\, \) = S(Cb — Cb) if A Cb = 0] to check 
the effect of uncertainties in the distribution of the bolometric corrections. 
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Note that inequalities (60) and (62) can be generalized to other bands (e.g., the X-ray band). Given the 
observed QSOLF and the bolometric correction in some other bands, it is worthy to make similar comparisons 
as is done in the B band below (see also discussion in § 5.2). 

4.2. Observational constraints on the nuclear/QSO luminosity evolution and BH growth 

In § 2.4, the model of the nuclear/QSO luminosity evolution is characterized by three parameters: the 
period of the nuclear activity in the first phase n, the characteristic luminosity declining timcscalc in the 
second phase td, and the efficiency e. In § 4.2.1, using inequality (60) or (62) (see also eq. 21), we compare 
Q SO (MB,to) with Tm b Joc3.i(Mb, to) obtained with different values of ri/rs p , td/ts p , and e to provide 
observational constraints on the QSO model and BH growth. In the figures of this section, the results are 
shown in the luminosity range —20 > Mb > —30. However, the results are fairly secure probably only in the 
range —23 > Mb > —26 (Boyle ct al. 2000) (which are marked by two vertical thin dot-dashed lines in the 
figures), and the results at Mb < —26 might be affected by the uncertainty of the small number statistics. 
In § 4.2.3, we show the effects of the uncertainty of the intrinsic scatter of the BH mass-velocity dispersion 
relation of local galaxies and the scatter of the bolometric correction of QSOs (Ai og M BH ,o an d ^c B )- 

4-2.1. ri/rsp, td/ts p , e, and the obscuration ratio 

In this subsection we always set Ai og M BH ,o = 0.27 dex (Tremaine et al. 2002) and Ac B — 4.3 (Elvis et al. 
1994). With the luminosity evolution described in § 2.4, below we study three models: in model (a), only 
the first (or "demand limited") phase is considered; in model (b), only the second (or supply limited) phase 
is considered; and in model (c), both the first and the second phases are considered. 

In model (a), td — and n is a free parameter (see eq. 30). We first show the result obtained by setting 
e = 0.1 in Figure 3 and then see the change of the result by changing the value of e in Figure 4. 

• In Figure 3a, Tm b ,iocai(MB, to) is shown as the dashed lines, and T^F* qso(Mb, is shown as the solid 
line; and their ratios 7\L opt (M,e, t ) are shown as the dashed lines in panel (b). As seen from Figure 3(a) 
and (b), T MB ^ oca \(M B , t ) and lZ opt (M B , to) increase with increasing ti/t Sp (= 0.3, 1, 4 from bottom to 
top at the faint end, respectively). The reason of the increasing is that increasing ti/ts p will increase 
the mass range (or the upper limit of the mass range) of local BHs that had nuclear luminosity L in 
their evolution history (see eqs. 18 and 36). In this model, we get the following constraints: 

(i) Constraints on t\/tq, p from the sensitivity of T/v/ B! i oca i (Mb, to) to ti/ts p : as seen from Figure 3(a), 
if Ti/rs p > 1, Tm b ,iocai(Ms, t ) becomes insensitive to the value of ti/ts p at Mb < —26. Our 
calculations also show that if ti/ts p > 4, 7m b , local (Mb, to) becomes insensitive to the values of 
ti/t Sp at M B < -23. For example, we have T Ms j ocal (n/r Sp = 10)/T MB;local (Ti/T S p = 1) - 1 < 
10% at M B < -26 and T MB;local (n/T Sp = 10)/T MB;local (ri/T Sp = 4) - 1 < 15% at M B < -23. 
The reason for this insensitivity is that the additional contribution to T MjJj i oca j (Mb, t ) at a given 
Mb due to the increase of ti/ts p comes from the BHs at the high-mass end (see eqs. 18 and 
36). Thus, for a given Mb, if t\/ts p is larger than a certain value, the masses of the BHs from 
which the additional contribution comes will be significantly high, and the additional contribution 
may be negligible since the local BHMF decreases sharply (or exponentially) at the high-mass 
end (Mbr\o >3x 10 8 M q ). The insensitivity of 7M B ,iocai(Ms, t ) to the value of ti/ts p above 
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suggests that if the real ti/ts p is larger than 4, it is more likely that only the lower limit of ti/ts p 
(~ 4 here), rather than an accurate value of ri/rs p , can be provided by fitting 7m b , local to) 
to TM B ,Qso(MB,to)- To get an accurate estimate of ti/ts p , it would require either observations 
at fainter luminosities or precise measurements of both the QSOLF and the local BHMF (e.g., 
with error much less than 10%) within the current luminosity range. The inscnsitivity of Tm b , local 
to ti/ts p still holds if td is non-zero (see eqs. 33 and 36). The suggestion above will also not be 
affected by changing the value of e, since as will be seen below (Fig. 4), changing e does not affect 
the shape of T MB ,ioca,\( M B, to) or TZ opt (M B , t ) for a given ri/r Sp . 

(ii) Constraints on ti/t Sp from the value of 1Z opt (M B , t ): Figure 3(b) shows that TZ opt (M B ,t ) is 
significantly smaller than 1 if Ti/rg p <C 1 (see the bottom dashed line with ti/ts p = 0.3); and 
TZ opt (M B ,to) can roughly exceed or be around 1 at —23 > M B > —26 if n/rgp > 1 (see the 
middle and top dashed lines obtained with ti/ts p = 1 and 4). Thus, to satisfy inequality (62), 
the QSO lifetime should be > ts p ~ 5 x 10 7 yr if e ~ 0.1. This conclusion will not be affected 
after setting a non-zero timescale td, since as will be seen below (Fig. 6), increasing td decreases 
TZ opt (M B ,t ) at least at M B 26. 

(iii) Constraints on the ratio of obscuration: as seen from Figure 3b, if ti/ts p > 4, lZ opt (MB, to) is 
about 2 at Mb — —23 and about 1 at Mb — —26. This result suggests that if ti/ts p > 4, there 
should exist some obscured QSOs (which are missed by optical surveys), and the ratio of obscured 
QSOs to unobscured (or optically bright) QSOs may be about 1 at the faint end (intrinsic absolute 
magnitude Mb ~ —23); compared to optically bright QSOs, few obscured QSOs are expected to 
exist at the luminous end (Mb ~ —26). If t\/t^ v ~ 1, few obscured QSOs/AGNs are expected 
to exist at both faint and bright ends. Note that these conclusions are obtained with e = 0.1 and 
td = 0, and may be affected if e > 0.1 or td > 0. 

• In Figure 4a and 4b, the dotted, dashed, and dot-dashed lines represent the results obtained with e = 
0.057 (the efficiency for the thin-disk accretion onto a Schwarzschild BH), 0.1, and 0.31 (the maximum 
efficiency allowed for the thin-disk accretion onto a Kerr BH with dimensionless spin parameter 0.998, 
see Thorne 1974; here we do not consider the possibility of extremely high efficiency, see, e.g., Narayan, 
Igumenshchev & Abramowicz 2003; Li & Paczyhski 2000), respectively. As in Figure 3, for each type 
of those lines, we show the results obtained with ti/ts p = 4, 1, and 0.3 (top, middle, and bottom lines, 
respectively). As seen from Figure 4, T MBt \ oca \(M B , f ) and lZ opt (M B ,t ) increase with increasing e 
without changing the shape (which can also be inferred from eqs. 31 and 35). We obtain the following 
constraints from Figure 4. 

(i) Constraints on e: as seen from Figure 4(b), if e = 0.057 (dotted lines), TZ opt (M B , to) are generally 
smaller than 1, which suggests that the efficiency of QSOs/AGNs cannot be significantly less than 
0.1. 

(ii) Constraints on ti/ts p : our calculations show that if e = 0.31, to satisfy inequality (62), it is 
required that ti/ts p > 0.2 or the QSO lifetime should be > 4 x 10 7 yr (see also the bottom 
dot-dashed line in Fig. 4b). 

(iii) Constraints on the ratio of obscuration: as seen from Figure 4(b), if e = 0.31 and ti/ts p > 4 

(see the top dot-dashed line), TZ opt (M B ,t ) is ~ 8 at M B 23 and ~ 4 at M B 26, which 

suggests that the upper limit of the ratio of obscured QSOs/AGNs to un-obscured QSOs/AGNs 
is - 7 at M B ~ -23 and - 3 at M B ~ -26. 
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The constraints obtained above from Figure 4 will not be affected by having a non-zero td, since as 
will be seen below (Fig. 6), increasing td decreases the Tm b ,\oc&\(Mb, to) or lZ opt (MB, to) at least at 
M B 26. 

In model (b), n = and td is a free parameter (sec eq. 32). Setting e = 0.1, the results obtained 
with td/tsp = 0.3, 1, and 4 are shown as dashed lines from bottom to top at the faint end, respectively, in 
Figure 5. As seen from Figure 5, for any given td, Tm b , local (Mb, to) is smaller than T^P* q SO (Mb, to) at 
Mb < —26 at least by a factor of ~ 10 if e = 0.1 (and at least by a factor of 2 or 3 if e = 0.31). This result 
suggests that to satisfy the expected relation between the QSOLF and local BHMF (inequalities 60 or 62), 
it is impossible to only have the luminosity-declining phase in the QSO luminosity evolution. 

In model (c), both t\ and td are not zero. Setting ti/ts p = 4, 1,0.3 and e = 0.1, we show the results 
obtained with td/ts p = 0.3 and 1 as the dashed and dotted lines, respectively, in Figure 6. The value of td 
affects the shape and value of Tm b ,\ochi(Mb, to) and TZ° vt (Ms,to)- As seen from Figures 3 and 6, increasing 
t d/t"Sp decreases Tm b , local (Mb, to) and lZ opt (MB,to) at high luminosities (e.g., Mb < —26). Figure 6b 
shows that at Mb ~ —26, lZ opt (MB,to) is significantly smaller than 1 (about 0.2-0.3 if td/ts p = 1 and 
about 0.3-0.6 if td/ts p = 0.3), which suggests that td should be significantly shorter than ts p , or the second 
(or supply limited) phase should not dominate the growth of local BHs (see eq. 33). Our calculation shows 
that this suggestion will not be significantly affected even if e = 0.31. 

By combining the constraints obtained in the three models above, below we summarize the results on 
the fundamental parameters of the QSO luminosity evolution and BH growth: 

1. The QSO mass-to-energy conversion efficiency is e > 0.1 (see Fig. 4). 

2. The period of the nuclear activity is longer than ~ rs p ~ 5 x 10 yr if e = 0.1 and longer than 
- 0.2r Sp = 4 x 10 7 yr if e = 0.31 (see Figs. 3 and 4). 

3. If the real n is larger than a certain value (~ 4rs p here), it is difficult to provide an accurate estimate 
on the value of t\ unless observations are extended to fainter luminosities or precise measurements of 
both the QSOLF and the local BHMF within the current luminosity range are available (e.g., with 
error much less than 10%; see Fig. 3). 

4. The characteristic declining timescale of the luminosity evolution in the second phase td should be 
significantly shorter than ts p (e.g., td < 0.3rg p if e = 0.1 and td < ts p if e — 0.31; see Fig. 6); thus, 
the second phase does not dominate the growth of local BHs (see eq. 33). 

5. There might exist a large number of obscured QSOs/AGNs, but the ratio of the obscured QSOs/AGNs 
to the unobscurcd (or optically bright) QSOs should be not larger than 7 at M B ~ —23 and 3 at 
Mb ~ —26 if e ~ 0.31 (see the top dotted line in Fig. 4) and not larger than 1 at Mb ~ —23 and 
negligible at M B 26 if e = 0.1 (see Fig. 3). 

According to the results above, the lower limit of the QSO lifetime (defined directly through the lumi- 
nosity evolution of individual QSOs) is ~ 4 x 10 7 yr. There are also other independent methods (e.g., Martini 
2003), such as the clustering of QSOs (Martini & Weinberg 2001; Haiman & Hui 2001) or the modeling of 
the QSOLF in hierarchical galaxy formation scenario (Kauffmann & Haehnelt 2000; Haiman & Loeb 1998), 
to provide constraints on the QSO lifetime (usually in the range of 10 6 — 10 8 yr). Here we do not make a 
comparison with those results since the meanings of the QSO lifetime in different methods or contexts are 
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Fig. 3. — Time integral of the observed LF of optically bright QSOs and the time integral of the QSOLF 
predicted from local BHs. The model of the luminosity evolution of the nuclear activity in § 2.4 is used. 
The curves are obtained by setting the parameters td = 0, e = 0.1, Ai og M B H,o = 0.27 dex, and Ac B = 4.3. 
(a) Dashed lines represent the time integral of the QSOLF predicted from local BHs T/Wb, local (Afe > *o) 
(eqs. 18, 34, 35 and 57); solid line represents the time integral of the observed LF of optically bright QSOs 
T MB,QSo( M B,t ) (eqs. 45 and 55). (b) Ratio of the time integrals n opt (M B ,to) (eq. 62). The horizontal 
thin dot-dashed line represents lZ opt (MB,to) = 1. The two vertical thin dot-dashed lines mark the region 
—23 < Mb < —26 in which the LFs of optically bright QSOs arc fairly secure (Boyle et al. 2000). In each 
panel, the dashed lines represent the results obtained with t\/tq, p = 0.3, 1, and 4 from bottom to top at the 
faint end, respectively. Inequality (62) implies that the dashed line should be consistent with the solid line 
in (a) or approximately equal to 1 in (b) if obscured QSOs are significantly less than optically bright QSOs, 
and the dashed line should be higher than the solid line in (a) or higher than 1 in (b) if obscured QSOs are 
significantly numerous compared to optically bright QSOs. To satisfy inequality (62), the QSO lifetime t\ is 
required to be > ts p ~ 5 x 10 7 yr (if e ~ 0.1). See more discussion in § 4.2.1. 
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Fig. 4. — Same as Figure 3, except that there are additional results obtained with e = 0.057 and 0.31 (not 
only e = 0.1 in Fig. 3) shown as the dotted and dot-dashed lines, respectively. The curves shift upward with 
increasing e. As seen from (b), generally or at least at Mb < —25, 7l opt (MB, to) obtained with e = 0.057 
(dotted lines) are smaller than 1, which suggests that e cannot be significantly smaller than 0.1. See more 
discussion in § 4.2.1. 
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Fig. 5. — Same as Figure 3, but with n = and td/ts p = 4, 1, and 0.3 (from top to bottom at the faint 
end, respectively; dashed lines; model [b]). (b) lZ opt (Mb ~ —26) is significantly smaller than 1 (which is 
still true even if e = 0.31). This figure suggests that to satisfy inequalities (60) and (62), it is impossible to 
only have the luminosity-declining phase (eq. 32, the second phase described in § 2.4) in the QSO luminosity 
evolution. See § 4.2.1. 



- 23 - 




Fig. 6. — Same as Figure 3, but with td/ts p = 0.3 (dashed lines) and 1 (dotted lines; model [c]). (b) 
TZ opt (Mb ~ —26) is significantly smaller than 1, which suggests that to satisfy inequalities (60) and (62), 
td should be significantly shorter than rs P (e.g., td < 0.3ts p ), or the luminosity-declining phase (the second 
phase described in § 2.4) should not dominate the growth of local BHs. See § 4.2.1. 



-24- 



not exactly the same. The upper limit of the lifetime can be usually constrained by the rising and falling of 
the characteristic luminosity of the entire QSO population as a function of cosmic time (roughly a few times 
10 9 yr; see eqs. 48 and 52). 

It is worthy to note that item 3 above would hold for many other methods to estimate the QSO lifetime 
by using the QSOLF, not only for the specific study in this paper, as a result of the sharp decrease of the 
QSOLF at the bright end and the limited luminosity range of the observations. 

Item 5 above does not exclude the possibility of the existence of a large number of obscured QSOs at 
high luminosities if the efficiency is high (e.g., e ~ 0.3), which is slightly different from the result obtained 
by comparing partial mass densities (see eq. 29) in Yu & Tremaine (2002) that obscured accretion is not 
important for the growth of high-mass BHs (> 1O 8 M ). The reason for the difference is partly because, in 
the present study, we adopt a high value of Ai og M B H o( = 0.27 dex) and include the effect of the scatter of the 
bolomctric corrections (for the detailed effects of Ai og M B Ho an d ^c B see § 4.2.3 and Figs. 8 and 9 below). 
Furthermore, Yu & Tremaine (2002) compare the partial mass density accreted in QSOs having luminosity 
higher than a certain value with the partial mass density in local BHs (see eq. 29), and the present study 
compares the number density (or the time-integral of the number density) at a given luminosity; hence, 
the result of the obscuration obtained in Yu & Tremaine (2002) is mainly for the ratio of all the BH mass 
accreted in obscured QSOs with intrinsic luminosity higher than a certain value, rather than the number 
ratio of obscured QSOs at a given luminosity in the present study. 

In addition, in Figures 3-6 (e < 0.31), inequalities (60) and (62) are always not satisfied at the bright 

end Mb < —29 26, which is mainly because the local BHMF declines more sharply than the QSOLF at 

the bright end (note that the velocity dispersion function of early-type galaxies is fitted by an exponential 
form at the bright end, but the QSOLF is fitted by a power law). The forms of the QSOLF and local BHMF 
at the bright end are probably affected by the uncertainty of the small number statistics; otherwise, the 

physical mechanism/parameters (e.g., QSO efficiency) of very luminous QSOs (Mb < —29 26) or the 

properties of nearby galaxies with very big BHs (> 10 9 M Q ) should be very different from those of the main 
population of QSOs or nearby early-type galaxies. 

4-2.2. Is it possible that most QSOs radiate at super- Eddington luminosities? 

In § 2.4 it is assumed that after the nuclear activity is triggered, QSOs first radiate at the Eddington 
luminosity for a period t\ and then radiate at sub-Eddington luminosities as a result of the decline of 
available accretion material supply. However, accretion with super-Eddington luminosities might occur in 
some cases. For example, when BH mass is small and the accretion rate is sufficiently high, the outflow 
pushed by the radiation may be trapped by the infalling gas and the energy is radiated away at a rate higher 
than the Eddington limit (see discussion in Blandford 2003), or when some strong density inhomogeneity is 
developed in a thin disk, the disk may also radiate at super-Eddington luminosity (Bcgelman 2002). 

Below we will see whether the expected relations given by equations (60) and (62) can be satisfied if 
all QSOs radiate at a luminosity higher than the Eddington luminosity, say, by a factor of I > 1. Note 
that in this case the value of the characteristic increasing timescale of the luminosity or BH mass, denoted 
by Tg p , is smaller than the value in equation (31) by a factor of I. In Figure 7 we show the time integral 
of the LFs obtained with n/r^ = 4, 1,0.3 (from top to bottom) and with I = 2 (dashed lines) and I = 5 
(dotted lines). The other parameters (such as td = 0, e — 0.1, etc.) in Figure 7 are the same as those in 
Figure 3. As seen from Figure 7(b), at -23 > M B > -26, either lZ opt (M B ,to) is smaller than 1, which is 
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inconsistent with inequality (62), or lZ opt (M b ,to) increases with increasing luminosity, which is inconsistent 
with current observations that few type II QSOs are observed and most observed obscured AGNs are at low 
luminosities (e.g., Barger et al. 2003). The results above suggest that it is unlikely that most QSOs radiate at 
a luminosity much higher than the Eddington luminosity (or accrete at an accretion rate much higher than 
the Eddington accretion rate) . We expect that accretion with supcr-Eddington luminosities operates maybe 
only at the very early stage of the nuclear activity, and will not significantly affect our results obtained in 
§ 4.2. 

4-2.3. Effects of the uncertainty of Ai og M BH , and Ac B 

The effect of the uncertainty of Ai og M B H,o is shown in Figure 8. In Figure 8 we show Tm b , local (Mb, to) and 
lZ opt (MB,to) obtained by setting Ai og M B H,o = (dot-dashed lines) and 0.4dex (dotted lines), as well as 
the results obtained in Figure 3 (shown by the dashed lines with Ai og M B H,o = 0.27 dex and the solid line). 
The other parameters of the dot-dashed and the dotted lines are the same as in Figure 3. As seen from 
Figure 8(b), the effect of the uncertainty of Ai og M BH o is negligible at Mb > —24; however, TZ opt (MB,to) 
or TM Bt iocai(MB, to) is significantly affected by the value of A} og m B h o a t Mb ^ —26. For example, at 
Mb = —26, lZ opt (MB,to) obtained with Ai og M BH , = 0.4 dex (dotted lines) is larger than that obtained with 
Ai og m B h.o — 0.27 dex (dashed lines) by a factor of ~ 2 and is larger than that obtained with Ai ogMBH = 
(dot-dashed lines) by a factor of more than <~ 4. 

The effect of the uncertainty of Ac B is shown in Figure 9. In Figure 9 we show the results obtained 
by setting Aq b — as the dotted lines, as well as those obtained in Figure 3 (dashed lines; Ac B = 4.3). 
The other parameters of the dotted lines are the same as in Figure 3. As seen from Figure 9(b), at the 
faint end (Ms < —24), lZ opt (MB,to) [or Tm b , local [Mb, to)] is insensitive to the value of Ac B ; however, at 
the luminous end (Mb < —26), lZ° pt (MB, to) obtained with Ac B — 4.3 (dashed lines) is larger than that 
obtained with Aq< b — (dotted lines; e.g., by a factor of about 3 at Mb = —27). 

The results above raise the importance of accurately measuring both Ai og m bh ,o an d ^C B and of studying 
the dependence of the bolometric correction on the BH mass, the accretion rate, or other physical parameters 
for precise understanding of the relation between the local BHMF and the QSOLF (especially at the luminous 
end, e.g., Mb < —26), the luminosity evolution of the nuclear activity, and the BH growth. 

5. Discussions 

5.1. The QSO luminosity evolution and the triggering history of the accretion onto seed 

BHs 

As part of the steps to understand the physics behind the QSO phenomenon and BH growth, it is important 
to investigate both the luminosity evolution and the triggering history of the accretion onto their seed BHs, 
which together control the shape and the value of the QSOLF (see eqs. 2 and 9). In the currently popular 
coevolution scenario of QSOs and galaxy spheroids in the cold dark matter (CDM) cosmology, it is generally 
assumed that QSOs are triggered by hierarchical (major) mergers of galaxies, and the triggering rate is 
controlled by the (major) merger rate of galaxies (or halos in less fine models; e.g., Kauffmann & Haehnclt 
2000; Haiman & Loeb 1998 and references therein). With an assumed luminosity evolution (e.g., usually a 
step function or an exponentially declining function in those models), the coevolution models can reproduce 
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Fig. 7. — Same as Figure 3, except that the QSOs are assumed to radiate at a luminosity higher than the 
Eddington luminosity by a factor of I (t Sp = Z -1 ts p , Ti/Tg p = 4, 1, 0.3, td = 0, and e = 0.1). The dashed lines 
show the results obtained with I = 2, and the dotted lines show the results obtained with 1 = 5. As seen from 
(b), at —23 > Mb > —26, lZ° pt (MB, to) cither is smaller than 1, which is inconsistent with inequality (62), 
or increases with increasing luminosity, which is inconsistent with current observations on obscured AGNs 
(e.g., Barger et al. 2003). This figure suggests that it is unlikely that most QSOs radiate at a luminosity 
much higher than the Eddington luminosity (or accrete at an accretion rate much higher than the Eddington 
accretion rate). See details in § 4.2.2. 
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Fig. 8. — Same as Figure 3, except that there are additional results obtained with Ai og M B H.o — an d 0.4 dex 
(not only Ai og M B H,o = 0.27 dex in Fig. 3) shown as the dot-dashed and dotted lines, respectively. The effect 
of the uncertainty of Ai og M B H,o i s negligible at the faint end (Mb > —24) but may be significant at the 
bright end (M B < -26). Sec § 4.2.3. 
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Fig. 9. — Same as Figure 3, except that there are additional results obtained with Ac B = (not only 
Ac B = 4.3 in Fig. 3) shown as the dotted lines. The effect of the uncertainty of Ac B is negligible at the 
faint end (M B > -24) but may be significant at the bright end (M B < -26). See § 4.2.3. 
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the observed QSOLF. However, in those models, as a result of many uncertainties in the estimate of the 
galaxy merger rate and the gas infalling rate for BH growth, it is hard to differentiate whether the change 
of the QSOLF is due to a change of the nuclear luminosity evolution or to a change of the triggering rate; 
hence, it is hard to give an accurate constraint on the QSO luminosity evolution. For example the QSOLF 
can be reproduced for a large range of the assumed lifetime (10 6 to 10 8 yr; e.g. Kauffmann & Hachnelt 2000; 
Haiman & Loeb 1998). 

In this paper, by investigating the relation between the QSOLF and the local BHMF, we separate the 
luminosity evolution of individual QSOs from the triggering rate of the QSO population, with the assump- 
tion that each local massive BH has experienced the QSO phases and BH mergers are ignored. As shown 
in equations (18) and (21), the triggering rate is circumvented, and only the luminosity evolution of indi- 
vidual QSOs is implicitly reflected by their lifetime 7iif c (MBH,o) and the luminosity probability distribution 
P(L|Mbh,o) hi their evolution history. Thus, the constraints on the QSO luminosity evolution (such as the 
QSO lifetime and the efficiency) obtained here do not depend on the poorly known QSO triggering history, 
and we expect that these obtained constraints could be further used to infer the triggering history of seed 
BHs or refine the coevolution model for QSOs and galaxy spheroids (see also Blandford 2003). 

5.2. Obscured QSOs/AGNs 

Based on the assumption that the extragalactic X-ray background is mainly contributed by QSOs/AGNs, 
a larger number of obscured QSOs/AGNs (more than the optically bright or unobscurcd QSOs by a factor 
of > 4) were expected to exist from the X-ray background synthesis model (e.g., Gilli, Salvati & Hasinger 
2001, and references therein). The existence of those sources is generally consistent with the expectation 
of the unification model of QSOs/AGNs (Antonucci 1993): if the dusty torus is along the line of sight, the 
light from these QSOs/AGNs will be absorbed by the dusty torus (this structure is proposed to explain the 
classification of Seyfert 1 and Seyfert 2 galaxies at low redshift z < 1) and these QSOs/AGNs might be 
missed from observations in the optical band. In this scenario, obscuration is generally a geometric effect, 
and the obscured fraction is determined by the opening angle of the dusty torus. If the opening angle is 
the same for all QSOs/AGNs, the fraction of obscured QSOs/AGNs will be independent of the intrinsic 
luminosity of QSOs/AGNs. 

Recent X-ray observations by Chandra and XMM-Newton confirm the existence of obscured AGNs/QSOs 
(e.g., Barger et al. 2003). However, the fraction of obscured QSOs/AGNs appears not constant, which is 
higher at low redshift (z < 1) than at high redshifts (z > 1.5), and the obscured QSOs/AGNs have smaller 
BHs compared to optically bright QSOs (e.g., less than a few times 10 8 M Q ; see a recent review by Fabian 
2003 and references therein). This observational result contradicts with the expectation from the simple 
unification model if the opening angles of the torus are the same for all QSOs/AGNs, which might suggest 
that either the opening angle of the torus is smaller in luminous QSOs/AGNs than in faint QSOs/AGNs or 
QSOs/AGNs are obscured only at the early stage of their nuclear activities (Fabian 1999). 

The relation between the local BHMF and the QSOLF established in this paper may provide a way to 
explore the underlying physics of the obscured QSOs/AGNs, independent of the X-ray background synthesis 
model. With the ongoing and future observations on obscured QSOs/AGNs (e.g., the X-ray deep surveys 
by XMM and Chandra), together with the observations on un-obscured QSOs/AGNs (e.g., by SDSS), we 
expect that more accurate constraints on QSO models, the BH growth, and the physical mechanism of the 
obscuration can be obtained, which will improve our understanding of the coevolution of QSOs and galaxy 
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spheroids. 

5.3. BH mergers after the quenching of the nuclear activity and BH ejections 

The effect of BH mergers on the local BH distribution function, which depends on the galaxy merger rate 
and the binary BH evolution process, is still very uncertain (e.g., Begelman, Blandford & Rees 1980; Yu 
2002; Milosavljevic & Merritt 2002; Volonteri, Haardt & Madau 2002). For simplicity, here we only discuss 
the effects of mergers after the quenching of the nuclear activity. Since mergers of small BHs form big BHs, 
BH mergers may make the BHMF increase at the high-mass end and decrease at the low-mass end. Thus, 
the time integral of the QSOLF predicted from the local BHMF (cq. 18) might be overestimated at the 
luminous end and underestimated at the faint end as a result of BH mergers. Indeed, the tentative result 
from the hierarchical galaxy formation scenario shown in Fig. 5 of Kauffmann & Hachnclt (2000) is that 
the BHMF at redshift is lower than the BHMFs at higher redshifts (0.5, 1, 2) by a factor of < 2 at BH 
mass < 10 8 - 5 - 9 Mq and is higher than those BHMFs at higher redshifts at BH mass > IO 85 - 9 M Q . Thus, at 
least the models ruled out by applying inequalities (60) and (62) at the luminous end (Mb ~ —26) without 
considering BH mergers will still be ruled out. 

Not all massive BHs may reside in galactic centers (see § 3.2.1 in Yu & Tremaine 2002). For example, 
BHs may be ejected from galactic centers through either interactions of three or more BHs (e.g., Valtonen 
1996) or gravitational radiation reaction during BH coalescence (Rees 2001), and BHs may also be left in 
galactic halos after galaxy mergers if the BH mass ratio of the two merging galaxies is small enough (e.g., 
< 0.001; Yu 2002). After considering the possibility that there might be some BHs that have experienced 
the QSO phases not locating in galactic centers, the time integral of the QSOLF predicted from local 
galaxies in this paper would increase. Currently, it is hard to give an estimate on the fraction of the BHs 
ejected from galactic centers or left in galactic halos during galaxy mergers by both theoretical models (also 
because both the BH and galaxy merger history and the BBH merger process are still very uncertain) and 
observations. Volonteri, Haardt & Madau (2002) study the assembly and merging history of massive BHs 
in the hierarchical models of galaxy formation and argue that the population of BHs wandering in galactic 
halos and the intergalactic medium at the present epoch contributes to the total BH mass density by < 10%, 
which (if true) suggests that ignoring BH ejections would not significantly affect the results of this paper. 

A complete and quantitative consideration of BH mergers or ejections is beyond the scope of this paper. 

6. Other possible applications 

In this section we discuss two possible applications of the work established in § 2, to the study of the demog- 
raphy of QSOs/AGNs and the demography of the hot stellar components of normal galaxies at intermediate 
redshift. 



6.1. Demography of QSOs/AGNs 

Consider such a galactic property V, which does not significantly change during the nuclear active phase and 
after the quenching of the phase and is closely correlated with the BH mass Mbh,o at present. Study of the 
relation between the BH mass and V in QSOs/AGNs, as well as in nearby galaxies, might provide valuable 
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information on QSO models and BH growth. Below we show a way to investigate the relation between the 
BH mass and V in QSOs/AGNs using the work in § 2. 

The posterior distribution of Mbh,o given the luminosity L of a QSO can be defined as follows: 

,o) n M BH (-^BH.O; tp) 

/ 00 dMBH,o/o dt/ t di i ^(i i ,M B H,o,i,t) ~ " T i;local (L,io) ' 

(64) 

where equations (17) and (18) are used. Thus, given the local BHMF and the luminosity evolution of 
individual QSOs £(Mbh,o> T )> which can be used to obtain riif e (MBH,o)-P(-^|-^BH,o), we may use equation 
(64) to get V(Mbh,o\L) and further use the local Mbh,o — V relation and get the distribution of V at a given 
L or Mbh in QSOs. Comparison of future observation results (see also current observational results on the 
BH mass-velocity dispersion relation in QSOs by Shields et al. 2003) with the expected relations between 
the galactic parameter V and the luminosity/BH mass of QSOs may further constrain the QSO luminosity 
evolution. 

In a separate paper (Yu & Lu 2004) we will explore the nuclear luminosity/BH mass and velocity 
dispersion relation in QSOs/AGNs in the way described above, assuming that the velocity dispersion of the 
hot stellar components of galaxies does not significantly change during the nuclear active phase and after 
the quenching of the phase. 



6.2. The distribution of velocity dispersions in elliptical galaxies and bulges of SO/spiral 

galaxies at intermediate redshift 

Study of the demography of galaxies at all redshifts may help us understand the formation and evolution 
of galaxies. The past decade has seen dramatic expansion of the knowledge on intermediate- and high- 
rcdshift galaxies in both observations and theories, as well as on local galaxies, such as the observations of 
the Lyman-break galaxies at redshift z ~ 2 — 4 (e.g., Adelberger et al. 2003), the semi-analytic hierarchical 
galaxy formation models to explain and predict the observational properties of both early-type and late- type 
galaxies (e.g., White & Frcnk 1991; Kauffmann, White & Guiderdoni 1993; Cole et al. 1994, 2000), etc. 
Below using the work in § 2, we point out a simple way to estimate the velocity dispersion distribution in 
elliptical galaxies and the bulge components of SO/spiral galaxies at intermediate redshift (1 < z < 2.5 here), 
which are still currently poorly known and are expensive and difficult to measure by observations. 

We define nM BH (^BH,o, t) as the BHMF of dead QSOs at cosmic time t so that nM BH (-^BH,0j £)dMBH,o 
represents the comoving number density of the BHs whose nuclear activities had been quenched before time 
t and the mass of these BHs is in the range Mbh,o — » Mbh,o + dMBH,o from the quenching time to the 
present time to. Similarly as the derivation to get equations (18) and (21), we have 

r-t f-oo 

/ * L (L,t)dt= / dM B H,or(M B H,o) J P(i|MBH ,o)^M BH (-^BH,o ) i), (65) 
Jo Jo 

as long as the condition 

t ft 
t> L {L,t)dt?s> / * L (L,t)dt (66) 

Jt-nite 

is satisfied. Condition (66) represents that most QSOs contributing to the integral of the QSOLF over the 
cosmic time — t have become quiescent (note the similarity to the necessary condition to hold eq. 18 in § 2 



L 
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that most or all QSOs are quiescent at present time to)- Using the observed QSOLFs shown in § 3.2, we 
find that condition (66) can be satisfied at redshift z < 2.5 if Tuf < 4rs p (e = 0.1) ~2x 10 8 yr. 

As shown in § 4.2, the comparison between relation (21) and observations may provide constraints 
on the QSO luminosity evolution or T(M BBfi )P(L\M BBfi ). Given t(M B h,o)P{L\Mbh,o) and the QSOLF, 
"m B h(-^bh,o, t) can be uniquely determined from equation (65) (which is just a typical inverse problem). 
According to Bayes's theorem, the distribution of the velocity dispersion of the hot stellar components of 
local galaxies given a BH mass M B h,o is given by: 

P(a\M BH ,o) = — 777 -t— • (67) 

We assume that the formation and the significant part of the evolution of bulges are simultaneous with 
the significant evolution (and/or formation) of their central BHs and that their velocity dispersions (and 
BH masses) do not significantly change after the quenching of the nuclei activity. Thus, there will be little 
evolution of the BH mass and velocity dispersion relation in those galaxies after they have experienced QSO 
phases, and the significant part of the evolution of the BH mass- velocity dispersion relation is recorded only 
in QSOs/AGNs. With P(a\M B H,o) obtained in equation (67), the velocity dispersion distribution of the hot 
stellar components of galaxies at low and intermediate redshift (z < 2.5) can be given by 

poo 

n a {a,t)= dM BHfi P(a\M BHfi )n MBH { M BH,o,t). (68) 
Jo 

Note that according to the assumptions above, for normal galaxies, the hot stellar components exist only 
in those galaxies containing massive BHs (BH ejections are ignored here), and those galaxies that have 
experienced QSO phases and contain BHs must have hot stellar components, which of course should be 
tested by future observations. The methods above would still be applicable even if not all of the bulges 
or galaxies containing BHs follow these assumptions, but as long as most of them do. In addition, if the 
formation of the hot stellar components occurs before the formation of central BHs or the QSO phase (that 
is, some hot stellar components may not contain BHs), the velocity dispersion distribution of the galaxies 
obtained by the method in this section will at least give the lower limit to their realistic distribution. 



7. Conclusions 

Assuming that each massive BH in nearby galactic centers has experienced the QSO phase and becomes 
quiescent at present, we have established a relation between the QSOLF and the local BHMF by studying 
the continuity equation for the BH mass and nuclear luminosity distribution and ignoring BH mergers. 
This relation compares the time integral of the QSOLF and that inferred from the local BHMF and only 
incorporates the luminosity evolution of individual QSOs. The triggering history of the accretion onto 
seed BHs is (implicitly) considered in the continuity equation but is circumvented in the relation between 
the QSOLF and the local BHMF. For comparison, the old relations between QSOs and local BHs on the 
total/partial BH mass densities (see eqs. 26 and 29; Soltan 1982; Yu & Tremaine 2002) include the effect of 
BH mergers, but the seed BH mass is ignored. The relation on the total BH mass density (eq. 26; Soltan 
1982) is unrelated with the luminosity evolution of individual QSOs, and the relation on the partial BH mass 
density (eq. 29; Yu & Tremaine 2002) assumes that the luminosity of QSOs is only an increasing function of 
their central BH mass (e.g., the Eddington luminosity in the calculation in Yu & Tremaine 2002). The new 
relation on the time integral of the QSOLF in this paper can be used to explore the luminosity evolution of 
individual QSOs (see § 1 and 2.3). 
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By applying observations into the relation established in this paper and assuming that the nuclear 
luminosity evolution includes two phases (first increasing at the Eddington luminosity with the BH growth 
and then declining), we find that the time integral of the QSOLF is generally consistent with that inferred 
from local BHs and obtain the following observational constraints on the QSO luminosity evolution and BH 
growth: (i) The QSO mass-to-energy efficiency e should be > 0.1 (see Fig. 4). (ii) The lifetime (defined 
directly through the luminosity evolution of individual QSOs) should be longer than ts p (~ 5 x 10 7 yr) if 
e = 0.1 and 0.2rs p (— 4 x 10 7 yr) if e = 0.31 (see Figs. 3 and 4). The characteristic declining timescale in 
the second phase should be significantly shorter than rs p , and BH growth should not be dominated by the 
second phase (when QSOs are accreting at sub-Eddington luminosities; see Fig. 6). (iii) The upper limit of 
the ratio of obscured QSOs/AGNs to optically bright QSOs is provided, which should be not larger than 7 at 
M B ~ -23 and 3 at M B ~ -26 if e = 0.31 and not larger than 1 at M B ~ -23 and negligible at M B ~ -26 
if e = 0.1 (see Figs. 3 and 4). (iv) It is unlikely that most QSOs are accreting at super-Eddington luminosities 
(see Fig. 7). The constraints above are obtained by assuming that the two accretion phases appear only 
once in the luminosity evolution, although this assumption is not required in the relation established in this 
paper. The possibility of more complicated accretion patterns deserves further investigation. 

We find that if the QSO lifetime is longer than a certain value (e.g., ~ 4rs p ; see Fig. 3), the time 
integral of the nuclear LF inferred from local BHs becomes insensitive to the value of the QSO lifetime, and 
thus it is difficult to provide an accurate estimate on the QSO lifetime unless observations extend to fainter 
luminosities or precise measurements of the QSOLF and local BHMF are available (e.g., with error much 
less than 10%). We also point out that this difficulty would also exist in many other methods to estimate 
the QSO lifetime by using the QSOLF, as a result of the sharp decrease of the QSOLF at the bright end 
and the limited luminosity range in observations. 

We show the importance of accurately measuring the intrinsic scatter in the relation between the BH 
mass and velocity dispersion of local galaxies and the scatter in the distribution of bolometric corrections 
of QSOs to precise understanding of the physics behind the QSO phenomenon and BH growth. Both of the 
scatters affect the shape and values of the time-integral of the nuclear LF, especially at the bright end. 

With the upcoming more precise measurement on QSOs (including both unobscurcd and obscured 
AGNs) and the demography of local BHs and galaxies (by SDSS, Chandra, XMM, etc.), the method pre- 
sented in this study would help to further explore the nuclear activity triggering and quenching mechanisms, 
obscuration of QSOs/AGNs, the demography of QSOs/AGNs, and the demography of normal galaxies at 
intermediate redshift and finally understand the physics behind the QSO phenomenon and the formation 
and evolution of galaxies. 

We thank Ravi Shcth for helpful communication on the velocity dispersion distribution of nearby galax- 
ies. We thank Norm Murray for helpful discussions. We thank the referee and Scott Tremaine for thoughtful 
comments. 
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